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abstract 


1 1 


The mechanism of transformation in shoaling of waves 
is one of the fundamental problems in coastal engineering. For 
efficient and reliable design of coastal structures it is very 
important to know the shallow water wave charac ter i s t i cs . In 
the present investigation an experimental study of the 
transformation of waves as they propagate from deep water to 
shallow water is presented. Regular waves and random waves are 
generated in a wave flume with a sloping beach in the 
laboratory. Random waves generated have JONSWAP spectra. 
Results show that the phase velocity measured is close to that 
obtained by using small amplitude wave theory. Essential 
differences with the linear dispersion relation have been 
found for random waves. The spectral peak for wave height 
decreases with decrease in water depth. This shows increasing 
nonlinearity with decreasing depth but is in agreement with 
published literature. The frequency corresponding to the 
spectral peak increases with decreasing water depth. 
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Chapter I 


INTRODUCTION 


Offshore technology has experienced a remarkable growth 
in application in the recent years. At present a wide variety 
of offshore structures are being used even under severe 
environmental conditions. While these are predominantly 
related to oil and gas recovery, they are also being used in 
other applications such as harbour engineering and ocean 
thermal energy conversion. The potential of a major 
catastrophic failure, both in terms of human safety as well as 
economic loss, underlines the critical importance of efficient 
and reliable design. 

The mechanism of transformation in shoaling of waves is 
a basic problem in coastal engineering. In fact, the action of 
waves is an important factor in the design of structures in 
the coastal and offshore zones. In traditional design 
iregular it ies of waves have not been taken into account 
directly. Irregular waves have been represented as a 
superposition of a few regular waves having a predetermined 
wave height and wave period. 



The present thesis is concerned with the study of 
shoaling of regular and random waves. For regular waves, 
change in the characteristics of waves such as wave height and 
wave length as the waves propagate on a sloping beach have 
been studied. For random waves the following questions have 
been addressed. 

<i) Change in the power spectra with water depth. 

<ii) Change in dispersion relation of spectral components in 
shallow water conditions. 

<iii) Change in the relaton between mean and peak 
frequencies of wave spectra propagating over sloping 
bottom. 

<iv) Variation of energy flux. 

In parallel with studies of Bendykowska and Werner<1988> 
on the transformation of shallow water wave spectra, the 
present thesis is an experimental study of transformation 
of shallow water wave spectra. 

The present thesis is organised along the 
following lines. Chapter II contains a review of the 
existing literature on wave shoaling. Chapter III describes 
experimental details including the equipments used, 

instrumentation and measurement. Results and discussion and 
conclusions are presented in Chapter V and VI, respectively. 



Appendix A presents the measured signals in graphical form. 
Appendix B contains the definitions of important functions and 
parameters and their method of calculation. Appendix C 
contains computer programmes for calculating these quantities. 



Chapter II 


LITERATURE REVIEW 


2.1 SHOALING OF REGULAR WAVES - As a wave train propagates 
into shallow water, its characteristics such as the wave 
height and wave length change. This process is described as 
wave shoaling. The shoaling effect may be estimated on the 
basis of a simplified linear wave theory. Commonly used 
assumptions for developing this theory are, 

<a) motion is two dimensional. 

<b) wave period remains constant. 

<c> energy flux in the direction of wave propagation is 
constant . 

These asumptions require that sea bed has 
gentle slope resulting in negligible wave reflection and that 
energy be neither supplied nor dissipated by wave breaking or 
bottom roughness. The assumptions hold good only upto the 
point of breaking. Shoaling curves derived using assumptions 
given above show the variation of H / Ho, C/ Co = <L/ Lo> and 
Cg/ C with h/ g T* given by linear wave theory, are given in 
Figure 2.1. Relationship similar to those given above can be 
derived on the basis of a finite amplitude wave theory, but it 
is doubtful if a single wave theory will apply over the range 





6 


of depths, that cover deep water and shallow water. 

Iversen (19?2> presented experimental data 

that shows the height of periodic waves on sloping bottoms to 

grow faster than that predicted by sinusoidal wave theory. 

Harihar Rama Ayyar<1970) has studied the shoaling of regular 

waves in water over a steep bottom. The factors affecting the 

breakers characteristics are the incident wave 

2 

steepness (H/gT ) and bottom slope <m). 

R. C. Nelson<1980) studied the effect of the bed 
slope on the characteristics of nonbreaking waves. The main 

results of his study are the following: 

<i) The characteristics of waves near breaking on a gentle 

slope cannot be same as those on a horizontal bed. 

<ii> The effect of bed slope of less than 0.1 on potential 

energy and wave celerity is minimal and for all practical 
purposes can be ignored. 

<iii) Bed slopes greater than or equal to 0.062 can 
significantly affect the ratio of crest height to wave height 
and trough depth to wave height. 

Ib A. svendson & 3 . Buhr Hansen (1976) have studied the 
periodic waves moving on a beach and have compared the results 
with the theories of sinusoidal and cnoidal wave shoaling. 

Their results can be summarised as follows: 

(i) Linear theory can predict shoaling as long as wave 

height to water depth ratio is small. 



<ij> At breaker point the wave length to water depth ratio 
appears to be independent of bottom slope for bottom slopes 
less than 1:10. 

(iii> Phase velocity followed the prediction of linear wave 
theory for h/ Lo > 0.1. Scattering may be present but can 

be attributed to small free secondary harmonic waves present in 
the f 1 ume . 

(iv) Wave height discrepancy is attributed to friction of 
the walls of the flume. 

2.2 MODELS FOR FREQUENCY SPECTRA - To study the random waves, 
mathematical models for wave spectra have been developed. 
These are in terms of a reference wind speed U appearing as a 
parameter. Of the following, Bretchneider and Pierson and 
Moskowitz spectra are most commonly used. JONSWAP spectra, an 
extension of the Pierson and Moskowitz spectrum produces a 
much sharper spectral peak, is more recent and involves 
additional parameters. 

Pierson— Moskowitz Spectrunif 1 964 ) -Th i s spectrum depends only on 
the wind speed and is given as.. 

V Z 7T / f 

Where a =8.1 x 10”® and is termed as Philip's constant. 

B is 7.4( g /2 „ u 3* 



JONSWAP Spectrum <1973>- It is the modification to the 
P 1 e r son-Moskow 1 tz spectrum to account for the fetch 
restrictions and to provide a sharply peaked spectrum. It is 
given as , 


S(f) = 


z 

g 


\- ^ r 


<2 n l"* 


L ^ ^ 

fo J J 

a = 

exp [ 

-< f-fo)^/ 

' 2 c" f" 

CX =: H 

a = 

o. 

.07 

for f S 


- 

.09 

for f 

peak frequency 

at which 

S<f) IS 


found to be related to fetch parameter as, 


fo = 2.84(gF*/l/) 

Parameters ao. and vh relate respectively to the width of the 
left and right side of specral peak. 

a is Philip’s constant and depends upon fetch parameter as, 
a = .066 (g F*/ l/]"' 


2.5 SHOALING OF RANDOM WAVES - 


transfer per unit 


width. 


The 

or 


average 

energy 


rate of 
flux. 


energy 
of a 



unidirectional random wave train is given as. 


P = (og Cg S(f) df 

Where Cg is group velocity. 

Thus transformation of a unidirectional spectrum 
during the wave shoaling is given by 

Cg S<f) = constant 

Testuo Sakai and Yuichi Iwagaki <1974) 
studied the transformation of random waves in shoaling water 
with the spectrum having one or two random peaks. It was found 
that the elementary peak power of random waves decreased with 
decrease in water depth. This reason can be explained 
qualitatively by theory of change of component wave height of 
finite amplitude waves in shoaling water. The damping due to 
internal viscosity and boundary friction was found to be 
negligible because it was of the same order of magnitude as 
the experimental error. 

Bendykowska and Werner<1988> in their investigation on 
the shallow water wave spectra studied the effect of 
nonlinearity on shallow water wave spectra. 

The waves were generated mechanically in a laboratory wave 
flume with a fixed bottom. Deviation from linear dispersion 
relationship was found, showing a vanishing dispersivity of 
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higher frequency spectral components. The mean frequency was 
seen to increase with decreasing depth. The relation of peak 
frequency to mean frequency varied in experiments from 0,9 to 
0,9 for deep and shallow water spectra respectively. 

The mam aim of present thesis is 
motivated by the work of Bendykowska and Werner < 1 988 ) , However 
it is an independent study of the behaviour of regular and 
random waves as they move towards the shore. 



Chapter III 


EXPERIMENTAL DETAILS 


5.1 EQUIPMEMT 


WAVEMAKER — The wavemaker is of the absorption type and 
comprises of following components. It has a flap which is 

wedge shaped paddle pivoted about its tapered end on a fabric 
hinge fixed to the bottom leading egde of the wavemaker box. 
The paddle is mounted in the box and extends the full width of 
the tank. The front of the box is sealed with a waterproof 
fabric gusset, behind which the paddle moves. This serves to 
prevent water entering the box with the result that the paddle 
operates only on the water in front of it. The paddle is 

driven by an electric servo-moter via a beryllium-copper drive 
belt which runs over a sector attached to the top of the 
paddle. The hydrostatic force on the flap is offset by 

springs. A piezo-electric transducer is sited between the 
sector and the paddle and senses the force acting on the 



paddle. Position is sensed by an incremental encoder mounted 
on the motor-shaft. 


CONTROL BOX — The control box contains the power amplifier for 
the wavemaker, along with the signal conditioning, force feed- 
back, wave absorption and nonlinearity correction circuity. 
There is also load protection circuit that disables the power 
amplifier under fault conditions. 

An encoder on the drive motor is used for 
position sensing. The signals from this are counted, then 
converted to give analogue position signal. This is filtered to 
produce the velocity signals required to optimise the paddle 
absorption characteristics. The digital signal from the counter 
chain is passed to an EPROM that is programmed to provide a 
square law function for hydrodynamic and hydrostatic second 
order correction. 


INTERFACE BOX — The interface box acts as the junction between 
the computer and outside world. All the signals to and from 
the interface card in the back of PC are run by a cable to the 
interface box. The inputs to this unit consist of 16 buffered 
analogue inputs with filters and four digital input circuits. 
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The analogue inputs can be independently viewed by switches in 
front panels. The outputs from the PC consist of two analogue 
signals for the wavemaker and four digital outputs. The 

wavemaking signals pass through a digital potentiometer to 
control the overall gain. Digial outputs are buffered by 

re lays . 

In addition, the interface box contains a sine wave 
generator for running the wavemaker without using the PC. The 
system is run by an internal power supply. The wavemaker 
control box return^ force .velocity and drive signals that 

are available on a 15 way D connecter on the rear of the 

interface box. 


INTERFACE CARD — The interface card is manufactured by 
Advantech as part PCL 81Z and plugs into an expansion slot of 
a PC. It provides two 12 bit analogue output channels, 16 
analogue inputs, 16 digital inputs and 16 digital outputs. 


WAVEGAUGE— This is a simple probefor sensing the elevation of 
water waves. It works on the principle of measuring the 
current flowing in a probe consisting of parallel stainless 
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steel wires. The probe is energised with a high frequency 
square wave voltage to avoid polarisation effects at the wire 
surface. The wire dips into water and current that flows 
between them is proportional to the depth of immersion. The 
current is sensed by an electronic circuit providing an output 
voltage proportional to the instantaneous depth of immersion 
and hence wave height. Itcan be used to drive a high speed 
chart recorder and/or a data logger. 

The standard form of the probe consists of a 
pair of stainless steel wire, 1.5 mm in diameter and spaced 
12.5 mm apart. However flexibilty of drive and sensing 
circuits enables a wide variety of probe configuration to be 
employed. 

The output voltage can be calibrated in terms of wave 
height by varying the depth of immersion of probe in still water 
by a measured amount and noting the change in the output 
signal. A special probe support member with accurately spaced 
positioning holes is available to facilitate this operation. 

The photographs of the equipments used in the 
present study are given in Photograph 1 to 5. 






Connection to tank experiment 




A / 


Photo. 1 



Photograph showing the wave flume 
with piping for filling and emptying 
the wave flume and wave monitor 



Photo, 4 


Photograph showing the storage 
oscilloscope and PC 


Photo. 5 



Photograph showing the wavemaker 
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5 . 2 SPECl FI CATI ONS- 

?.2.1 OSCILLOSCOPE 
Model number 

Manufacture 

Number of input channels 
I nput coupl 1 ng 
Storage mode 

Normal mode 

Vertical sensitivity 

Vertical gain accuracy <dc> 
Time base accuracy 
Time base range 
Maximum sample rate 
Memory depth 

Signal averaging 

3.2.2 PC 
Type 

Manufacturer 


3.2.3 INTERFACE CARD 
Mode 1 


: HP 34301 A 

Digitizing oscilloscope 
: Hewlwett Packard 
: 4 

: AC . DC 

: Oscilloscope digitizes the 
input signal before 
dispaying the trace 
: Behaves like a conventional 
analog instrument 
: Maximum 3 mv/div 
Minimum 3 v/div 
: + 1.5 « 

: 0.005 * 

: 2 ns/div to 3 s/div 
: 10 Msa/s 

: 301 Points (display) 

1024 Points (via HP-IB) 

: 1 to 1024 


: 386 SX 

! WIPRO (PVT) Ltd. 
Banga lore 


: PCL 812 
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5.2.4 WAVE GAUGE AND MONITOR 
Type 

Manufacturer 
Output voltage 
Range 

Frequency response 

5.2.5 WAVEMAKER 
Type 

Manufacturer 

Component parts 


Power supply 


: Cap 1 c 1 tance 

: Armfield Ltd. U.K. 

: + 10 volt 

: 500 mm 

: 10 Hz 

: Absorption wavemaker 

: CAN Electrical Industries 

Ltd, England. 

: (a) Paddle assembly 
(b> Hinge 

(c) Membrane seal 

(d) Coupl i ng wire 

(e) Spring 

(f) Anchor chord 

(g) Motor tacho 

(h> P.C.B. (Contrl board) 

: +7.5 V and -7.5 V, 250 ma and 
+15 V and -15 V, 10 Amp per 
paddle system 


! 
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j . ? prcx:edure 

Experiments were conducted in a wave flume which 
is 24 meters long. 0.45 meter wide and 0.90 meter deep. A wave 
generator of flap-hinge type was installed at one end of the 
flume. A uniform sloping beach was set in the flume. This 
beach was made out of aluminium plates. The slope of the beach 
was 1 / 20. A schematic sketch of the wave flume is given in 
Figure 5.1 . In all experiments water depth in the constant 

depth region was kept at 70 cm. Positions ofPC, interface box 
and control box were connected as given in the Figure 5.2. 

First of all, wavegauge was calibrated by 
raising or lowering it by a known amount in still water. This 
operation was facilitated by mean of a calibrated holder which 
has series of holes drilled along its length, accurately 
spaced at 10 mm. The probe was raised or lowered by a multiple 
of 10 mm by using the pin provided as a stop. 

After calibrating the probe, regular and 
random waves with a JONSWAP spectrum were generated by PC 
using the computer programmes given in Appendix C. For 
regular waves, data was collected at 10 places in constant 
water depth and at 80 cm interval on the slope upto a depth of 
5 cm. For random waves, experiments were conducted for three 


different cases. These cases are described in Table 5.1. Data 



was collected at four locations in constant water depth and at 
1 meter interval on the slope upto a depth of 5 cm. 

Regular waves were generated with 1 Hz frequency and 5 
cm amplitude. Three JONSWAP spectra waves were generated with 
a peak frequency of 1 Hz, 0.75 Hz and 0.60 Hz and peak 
enhancement factor of 6.0, 1.0, and 1.0 respectively. 

The signals were collected by the PC in terms 
of voltages. These were converted to water surface 
elevation by using the calibration curve. 

Spectral analysis of the water surface 
elevation signal was performad by using FFT software on HP 
main-frame computer. The phase difference 4 > between signals of 
the two probes was obtained from cross spectral analysis. This 
phase difference was used to determine the wave number and 
celerity. Data having the cross correlation coefficient factor 
of greater than 0.8 has been identified as acceptible. Group 
velocity and energy flux for regular waves were calculated by 
using small amplitude wave theory. 



Where 


2A 


Table 3.1 parameters of JONSWAP SPECTRA 


S(f) = 


a g 




5 r f -v-* ■ 
exp|^- _ ) 


a = exp 


(f-foY/ 2<r*fl 


] 


O' =0.07 

O. 

for 

f < 

f O 

= 0.09 

V 

for 

f > 

f o 


g = 9.81 m/s^, a = Philip's constant = 0.0081 



f O 

CHs) 

r 

CASE 1 

1.00 

6.0 

CASE 2 

0.75 

1.0 

CASE 3 

0.60 

1.0 
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Chapter IV 


RESULTS AND DISCUSSION 


Results have been presented seperately for regular and 
random waves . 

REGULAR WAVES- 

Wave height and root mean square wave heght have been 
tabulated in Table 4.1. Figure 4.1 shows the recorded 
variation of wave height as a function of the local depth 
parameter h/Lo for a deep water steepness Ho/ Lo =0.055. It 
can be observed that the wave height decreases as the depth of 
water decreases. The experimental data follows the small 
amplitude wave theory for h/Lo larger than 0.15. For h/Lo < 
0.15, experimental data does not follow the small amplitude 
wave theory. 

The dimensionless phase velocity. group velocity 
normalised by phase velocity and energy flux are given in 
Table 4.2. The variation of dimensionless phase velocity and 
Cg/C with the dimensionless depth h/Lo are shown in Figure 4.2 
and Figure 4.5 respectively. The data is seen to follow the 
trend predicted by small amplitude wave theory. The variation 
of Cg/C with h/Lo shows that this increases as the water depth 
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decreases. Figure 4.5 shows the variation of the energy flux 
with water depth. Energy flux is seen to decrease with 
decreasing depth. This is against the assumption of a constant 
energy flux used in small amplitude theory. 

The deviation of experimental results from 
theoretical predictions can be attributed to the following 
factor . 

Propagation of energy flux to the water below the beach. It 
is believed that this factor is the most significant one. 

RANDOM WAVES- 

Cross-spectral analysis was applied to two 
simultaneous records taken 10 cm apart. Wave numbers K(<o> and 
celerity C(w) of the spectral components for each frequency 
were determined. The effect of nonlinearity was examined by 
varying the peak frequency of JONSWAP spectra in the range of 
5.76 < wp < 6.28 (rad/sec) at different water depths. 

Power spectra of wave heights measured 
at different depths are presented in Figure 4.6 to Figure 4.8. 
There is a substantial drop in spectral peak . This shows the 
increasing nonlinearity of wave motion at small depths. Figure 
4.6 to Figure 4.8 explain the increasing nonlinearity of wave 
motion on small depth. 

Figure 4.9 to Figure 4.12 show a comprehensive plot of 



wave height spectra for the three cases summarised in Table 
5.1. Results have been shown in constant depth region followed 
by the spectra for variable depth portion of the flume. The 
plots are practically identical in constant depth region. As 
one moves over the slope, wave having a given value of y 
develop an increase in spectral peak with decrease in peak 
frequency . 

From Figure 4.1> to Figure 4.16 results of wave 
frequency u> against the wave number K with the theoretical 
curve are plotted. All measured results indicate that the 
pridiction of linear theory are true only in the vicinity of 
peak frequency. but as the depth of water decreases 
nonlinearity also increases because we get more deviation from 
the theoretical curve. 

From Figure 4.17 to Figure 4.19 the line C({o)/Co(w) = l is 
the theoretical relation. The phase velocity in range ojp < w < 
2ojp follows the theoretical curve. Most of spectral components 
of to > Ztop travel at higher velocity than the theoretical. The 


presence 

of 

bounded 

wave 

1 s one 

of 

the 

expl anat i on . 

Same 

results 

are 

obta i ned 

by 

Crawford 

et 

al 

in 1981 in 

the i r 


numerical study. 

Results concerning the changes of mean frequency of 
the spectrum, when propagating from deep water (too) over 
decreasing depth to are shown in Figure 4.20. This figure shows 
an increase of mean frequency when water depth decreases. 
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Changes of relationship between the peak top and the mean 
frequency to of the spectrum as the function of relative water 
depth Kpo . h <Kpo = tn^/g) is given in Figure 4.21. The relation 
top/ to varies from 0.5 on deep water to 0.2 for very small 
relative depth. These changes are due to change of the shape 
of the spectra. 

Table 4.5 to Table 4.5 show energy flux, spectrum 
width parameter and coherence factor. The coherence factor is 
quite high (greater than 0.8) in ma)ority of experiments. 
Spectral width parameter is about 0.9 in all experiments, 
showing that generated random waves are practically Gaussian. 

The energy flux calculations employ the experimental wave 
height spectra. Wave number calculated on the basis of the 
dispersion relation produce the column of energy flux marked 
Theor. Energy flux corresponding to peak frequency is marked 
P(top> and that averaged over the range 0 S o>p > 5 is markec 
P<ave). The experimental energy flux calculation employes 
experimentally determined values of the wave number. These 
calculations show the energy flux to behave erratically and 
has not been used for analysis. The theoretical values 
<P<Theor>> show a reduction in energy flux over the beach. 
This is most prominent for depths less than 57.5 cms . As in 
the case of regular waves, this is attributed to the 
propagation of the energy flux to the water below the beach. 



Table 4.1 calculation of wave height and r. m. s. wave 
height for regular waves 


Hoi 

= 8.34 cm , 

Hc.2 = 9. IE 

tern f = 

1 .0 Hz 



s. 

No. 

distance 
from wave- 
maker 
( m ) 

6 /La 

Hi/ Hoi 

Hz/ Hoz 

n^/ Hoi 

j 

n^/ Hoz 

1 

4.0 

0. 4487 

0.9846 

0.9846 

0.349 

0.337 

1 

7 

5.0 

0.4326 

1.0190 

0.9836 

0.340 

0.334 

j 

3 

5.8 

0. 4075 

1.0770 

1 . 0570 

0.365 

0.356 

4 

6.6 

0. 3814 

0.9832 

0.9771 

0.340 

0.337 

5 

7.4 

0.3557 

0.9520 

0. 9705 

0.331 

0.342 

6 

8.2 

0.3301 

0.9580 

0. 9564 

0.330 

0.328 

7 

9.0 

0. 3044 

0.9221 

0.9117 

0,316 

0.313 

8 

9.8 

0.2788 

0.9410 

0.9204 

0.316 

0.313 

9 

10.6 

0.2532 

0.9377 

0.8856 

0.323 

0.305 

10 

11.4 

0. 2275 

0.9017 

0.8801 

0.310 

0.301 

11 

^ /-k /*n 

X X. ■ .£1 

0. 2019 

0.8934 

0.8447 

0.351 

0.288 

12 

13.0 

0. 1762 

0.8742 

0.8137 

0.286 

0.274 

13 

13.8 

0. 1506 

0.8038 

0.7581 

0.286 

0.264 

14 

14.2 

0. 1250 

0.7775 

0.8000 

0.251 

0.263 

15 

15.0 

0.0993 

0.7341 

0.7603 

0.255 

0.260 

16 

15.4 

0. 0865 

0.7440 

0.7734 

0.242 

0.268 

17 

15.8 

0.0737 

0.7544 

0.7309 

0. 228 

0.231 

18 

16.2 

0.0608 

0.6814 

0.5577 

0.215 

0.202 

19 

16.6 

0.0480 

0.5389 

0.4117 

0. 195 

0. 167 

20 

17.0 

0.0352 

0.6241 

0.5034 

0. 163 

0. 134 
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Table 4,2 Calculation of C/Co . C/Cg, Pi and Pz 
for regular waves 


Co-Lo = 1.56 meter f= i . o Hz 


m 

distance 
from wave - 
maker 
(m) 

a 

> 

0 

C/tg 

n 

0 

P 1 

(N/ Sec) 

(N/ Sec) i 

1 

4.0 

0.4487 

0.511 

0. 8910 

5,90 

7. 12 


5.0 

0.4326 

0.513 

0. 8884 

6.27 

7.06 


5.8 

0.4075 

0.516 

t 

0.8717 

6.96 

8. 10 


6.6 

0.3814 

0.517 

0.8269 

5.50 

1 

6.56 


7.4 

0.3557 

0.526 

0.8397 

5.33 

6.70 


8.2 

0.3301 

0.532 

0. 8205 

5.33 

6.42 

7 

9.0 

0.3044 

0.538 

0. 7884 

4.80 

6.24 

8 

9.8 

0.2788 

0.559 

0.8141 

4.99 

6. 19 


10.6 

0.2532 

0.584 

0.8333 

5.70 

6. 14 


11.4 

0.2275 

0.569 

0.6987 

4.30 

4. 96 


13.0 

0. 1762 

0.655 

0.7500 

5.00 

5.73 


13.8 

0. 1506 

0.634 

0.6000 

3.27 

3.51 

13 

14.6 

0. 1250 

0.735 

0.6858 

4.05 

5. 18 

14 

15.0 

0.0993 

0.767 

0.6070 

3.33 

4.30 

15 

15.4 

0.0890 

0.795 

0.5814 

3.31 

4.27 

16 

15.8 

0.0737 

0.868 

0.6600 

4.21 

4.96 
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Table 4.3 Calculation of coherence foctor, spectrum 
width parameter and energy flux for random 

waves < CASE 1 > 


s . 

No. 

distance 
from the 
wavemaker 

< Tn> 

depth 

of . 
water 

< cm > 

coher - 

ence 

factor 


Pttop) 

N/ sec . 
(Theov) 

ptave) 

N/sec . 
(TheoV) 

P(.0>p) 

N/sec. 
(Exp ) 

P(ave) 

N/sec. 
(Exp ) 

1 

2 . 5 

70.0 

0.80 

0.90 

2.60 

1.13 

0.05 

0.513 

2 

3.0 

70.0 

0.81 

0.90 

2.64 

1.13 

0 . 15 

0.048 

3 

3 . 5 

70.0 

0.79 

0.90 

2.65 

1.12 

2.24 

0.994 

4 

4 . 0 

70.0 

0.84 

0.90 

2.67 

1.12 

1 

0 . 15 

0.085 

5 

5 . 0 

67.5 

0.81 

0.88 

2.53 

1.07 

2.43 

1.410 

6 

6.0 

62.5 

0.81 

0.90 

2 . 53 

1.07 

0.08 

0.401 

7 

7.0 

57.5 

0.83 

0.88 

o Q«=: 

1.01 

0.08 I 

0.071 

8 

8.0 

52.5 

0.81 

0.90 

2 . 25 

0.94 

2.45 

1.060 

9 

9.0 

47.5 

0.80 

0.88 

2 . 12 

0.87 

0 . 15 

0.084 

10 

10.0 

42.5 

0.82 

0.88 

1.87 

0.80 

2.67 

0.989 

11 

11.0 

37.5 

0.84 

0.88 

1.65 

0.72 

2.91 

1 . 110 

12 

12.0 

32.5 

0.79 

0.87 

1.58 

0.64 

0.70 

0 . 163 

13 

13.0 

27.5 

0.81 

0.87 

1.30 

0.54 

2.60 

1.030 

14 

14.0 

22.5 

0.72 

0.87 

1.07 

0.44 

0.07 

0.075 

15 

15. 5 

15.0 

0.67 

0.87 

0.62 

0.30 

0 . 12 

0.257 

16 

16.0 

12.5 

0.98 

0.88 

0.59 

0.25 

1.67 

0 . 690 

17 

16.5 

10.0 

0.98 

0.90 

0.46 

0.20 

1.80 

0.794 
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Table 4.4 Calculation of coherence foctor, spectrum 

width parameter and energy flux for random 

waves (CASE 2) 


s . 

No. 

d ist ance 
from the 
wavemaker 

< Tn> 

depth 

water 

< CTn> 

coher - i 
ence 

factor 

e 

PtWp) 
N/sec . 
(Theor ) 

PCave) 

N/sec . 
(Theor ) 

P(«p) 

N/ sec . 
(Exp ) 

P(ave) 

N/sec . 
(Exp ) 

1 

2 . 5 

70.0 

0.85 

0.92 

2 . 19 

1.86 

0.09 

0 . 185 


3.0 

70.0 

0.86 

0.92 

2 . 33 

1.91 

0.08 

0.204 

3 

3 . 5 

70.0 

0.85 

0.92 

2 . 18 

1.89 

1.92 

1.57 

4 

4.0 

70.0 

0.88 

0.92 

2.33 

1.90 

0.07 

. 205 

5 

5. 0 i 

67.5 

0.86 

0.91 

2.24 

1.83 

2 . 19 

1.70 

6 

6.0 

62.5 

0.86 

0.92 

2 . 13 

1.83 

2.00 

1.60 

7 

7 . 0 

57.5 

0.87 

0.91 

2.01 

1.76 

0.09 

0.08 

8 

8.0 

52.5 

0.85 

0.91 

1.92 

1.61 

0.09 

0 . 131 

9 

9 . 0 

47.5 

0.85 

0.91 

1. 89 

1.51 

0.06 

0.233 

10 

10.0 

42.5 

0.86 

0.91 

1.63 

1.35 

1.96 

1.660 

1 1 

11.0 

37.5 

0.86 

0.91 

1.41 

1.19 

2 . 03 

1.780 

12 

12.0 

32.5 

0.83 

0.91 

1 . 19 

1.06 

0 . 10 

0 . 183 

13 

13 . 0 

27.5 

0.84 

0.90 

1.00 

0.879 

1.84 

1.500 

14 

14.0 

22.5 

0.76 

0.90 

0 . 81 

0.703 

0.03 

0 . 126 

15 

15.0 

17.5 

0.74 

0.90 

0.58 

0.554 

0.06 

0 . 197 

16 

15.5 

15.0 

0.73 

0.90 

0.58 

0.474 

1.26 

1.060 

17 

16.0 

12.5 

0.98 

0.91 

0 . 44 

0.389 

0.07 

0.370 

18 

16.5 

10.0 

0.98 

0.91 

0.35 ! 

0.318 

0.07 

0.385 















able 4.5 Calculation of coherence foctor, spectrum 

width parameter and energy flux for random 

waves (CASE 5) 


s. 

No. 

distance 
from the 
wavemaker 

< m > 

depth 

of, 

water 

< c m > 

coher - 

ence 

f ac t or 


PtC^p) 

N/sec . 
(TheoV.) 

ptave) 

N/sec . 

( Theoy.> 

p(.wp) 

N/sec . 
(Exp ) 

ptave) 

N/sec . 
(Exp ) 

1 

2.5 

70.0 

0.91 

0.94 

4.20 

3.040 

4.78 

2.660 

■r> 

3.0 

70.0 

0.91 

0.94 

4.60 

3.060 

0.04 

0.214 

3 

3.5 

70.0 

0.79 

0.93 

1.98 

2.030 

1.09 

1.830 

4 

4.0 

70.0 

0.93 

0.94 

4. 16 

3.080 

0.09 

0.311 

5 

5.0 

67.5 

0.92 

0.94 

4.57 

2.980 

5.24 

3.721 

6 

6.0 

62.5 

0.91 

0.94 

4.02 

2.950 

3.72 

4.962 

7 

7.0 

57.5 

0.92 

0.94 

4.52 

2.870 

5.83 

2.640 

8 

8.0 

52.5 

0.91 

0.93 

3.42 

2.610 

3.06 

2.546 

9 

9.0 

47.5 

0.91 

0.93 

3.47 

2.450 

0.03 

0.356 

10 

10.0 

42.5 

0.91 

0.93 

3.02 

2.200 

0.07 

0.033 

11 

11.0 

37.5 

0.90 

0.93 

2.43 

1.910 

0.08 

0.387 

12 

12.0 

32.5 

0.88 

0.93 

2.31 

1.710 

0.03 

0.059 

13 

13.0 

27.5 1 

0.89 

0.93 

1.88 

1.460 

2.57 

3.660 


14.0 

22.5 

0.82 

0.93 

1. 49 

1.150 

1.83 

1.540 


15.0 

17.5 

0.82 

0.93 

1. 10 

0.890 

1.52 

1.530 

16 

15.5 

15.0 

0.79 

0.93 

0.69 ' 

0.772 

1.76 

2.322 

17 

16.0 

12.5 

0.97 

0.94 

0.68 

0.606 ! 

0.03 

0.328 

18 

16.5 

10.0 

0.97 

0.94 

0. 29 

0.450 

0. 18 

0.298 
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Chapter V 

CONCLUSIONS 


The following conclusions have been drawn in the present 
study . 

<1> For regular waves the experimental results regarding 
phase velocity for slope 1:20 are quite close to that of 
theoretical results by small amplitude wave theory. There is a 
considerable deviation in terms of wave height. 

(2) <i> For random waves a departure from the theoretical 

dispersion relation is found in the experiments. This 
departure increases with decrease in water depth. 

<ii) The mean frequency w increases with decreasing water 

depth . 

<iii) The peak of power spectra of wave height decreases 


as water depth decreases. 
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APPENDIX B 

DEFINITIONS OF IMPORTANT PARAMETERS AND METHOD OF THEIR 

CALCULATION 


WAVE CELERITY— is velocity of individual harmonic waves. It is 
calculated by using the phase difference betwwen the two 
signals reaching the two probes at any time. It is denoted by 
C. 

WAVE NUMBER is number of waves present per unit length of 
wave. It is given as, 

K = w/ C 

GROUP VELOCITY is the speed of propagation of afinite group 
of waves. This is the velocity by which energy is transferred 
in the direction of wave propagation. When calculated using 
small amplitude wave theory, one gets, 

Cg = 0.5 ^ 1 + <2Kh>/ sinh<2Kh>jc 

ENERGY FLUX is the average rate of energy transfer per unit 
Width acoss a plane of constant cross-section. 



6J 


For regular waves it can be calculated by using small 
amplitude wave theory from the formula. 


P = E.Cg 


Where E = Average energy density 
= P g 


For random waves 


P = p g Cg S<f) df 


SPECTRAL WIDTH PARAMETER is a measure of non-randomness of a 
irregular wave and is given by 


Where 


m i/2 



O 4 


m 

n 


S(f) df 


PROBABILITY DENSITY FUNCTION - For any stationary signal 


x(t> 
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P<x) = 

Li m 

T CO 

[ Tx/ < A X T 

>] 



where Tx is 

the 

time during 

wh 1 ch 

the 

s i gna 1 

amplitude lies between 

X and 

(x +Ax > and 

T 1 s 

tota 1 

s i gnal 


time. 

For a random signal P<x) tends towards some limited 
value as T approaches infinity. 

P(x> = r ( Tx / Ax 


Power spectral density function is a measure of energy 

contained within the signal in each frequency. 

If x( t.f, Af > IS that part of signal x(t> which lies in 
frequency band Af centred at frequency f. then the mean square 
value within this band is 

S*(f) = f-i- x’*' X* <t. f. Af> df ] 

T -> 00 t ^ o J 

The power spectral density can be defined by 

E<f> = [ S^(f> / Af X* 1 

Where x* is the total mean square. In practice, E<f> is 
found using Fourier transform. 
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APPENDIX C 

LISTING OF PROGRAMMES 


CONVERSION OF VOLTAGE SIGNALS TO ELEVATION SIGNALS 

dimension al ( 11 00 > , a2 ( 1 1 00 > , bl ( 1 1 00 ) . b2 ( 1 1 00 ) . 1 1 me ( 1 1 00 ) 
common / calib / v ( 40 . 2 ) . h< 40 , 2 ) . nl , n2 

print * , 'specify number of calibration points for probes 1 & 2 ' 
read * , nl , n2 

open(unit=12.file='calib.dat*,status='old’> 
read< 12 .*> (v(i.l>,h(i.l).i=l.nl) 
read( 12 .*) (v(i.2),h<i.2).i=l.n2) 
close<unit=12,status=’keep’ ) 

nu = l 0 

n=2**nu 

pi =3. 141 3926 

open<unit=8,fi le=’p.dat' ,status=*old' ) 
read(8.*> ( t ime < i ) , al < i ) , i =1 ,n) 
read(8,*) (t ime < i ) ,a2< i ) , i=l ,n) 
c 1 ose ( un i t=8 , status= ’ keep ' > 

do 10 i = 1 , n 
al<i>=10.0*al<i) 
a2<i)=10.0*a2<i ) 
bl < i)=convert<al<i ) ,1) 
b2< i )=convert<a2< i ) ,2) 
cont i nue 

open(unit=16,fi le=’ph.dat' ,status=' unknown ' > 

write(16.111) <time(i>.bl(i),i=l,n) 

format <lx,fl2.3.8x,fl2.3) 

wr i te( 16. 211 > ( t ime < i ) . b2 ( i ) . i =1 .n) 

format <lx,fl2.3.8x,fl2.3> 

close(unit=16 . status = ' keep ’ > 

stop 

end 

function convert <x , i p) 

common / calib / v< 4u , 2 ) ,h(40 ,2) ,nl ,n2 

i f < i p . eq . 1 ) n=nl 
i f < i p . eq . 2 ) n = n2 
dol0i=l,n-l 

i f ( X . 1 1 . V ( i , 1 p) . and . X . ge . V ( i + 1 , i p) > go to 100 
cont i nue 

print *, 'interpolation failure' 
print 'probe number', ip 

print X 

gra§=<h( i+l,ip)-h<i.ip>)/<y<i+l,ip)-v(i.ip)) 
convert = h< j , ip> + grad*<x-v( i , ip)> 


return 

end 



APPENDIX C 


LISTING OF PROGRAMMES 


C CONVERSION OF VOLTAGE SIGNALS TO ELEVATION SIGNALS 

dimension al ( 11 00 ) . a2 < 1 1 00 ) , bl < 1 1 00 > . b2 < 1 1 00 ) . t ime ( 1 1 00 > 
common / calib / v < 40 . 2 ) . h ( 40 , 2 > , nl , n2 

print *, 'specify number of calibration points for probes 1 & 
read * , nl , n2 

open<unit=12,fi le=’cal ib.dat' ,status='old' ) 
read<12 .*> < v < i . 1 ) . h ( i . 1 ) . i =1 . nl ) 
read<12 .*> <v(i.2).h<i.2).i=l.n2) 
close(unit=12,status=’keep' ) 

nu= 1 0 

n=2**nu 

pi = 3 .1415926 

open<unit=8,fi le= 'p.dat' ,status= 'old' ) 
read<8,*> < t i me < i ) , al < i ) . i =1 , n) 
read<8,*) < t i me ( i ) . a2 < i ) , i =1 , n> 
close(unit = 8.status = 'keep' > 

do 10 i =1 . n 

al< 1 )=10.0*al < i> 

a2< i )=10.0*a2( i > 

bl < i ) = convert (al ( i > , 1 > 

b2< i )=convert<a2< i ) .2) 

10 continue 

open<unit=16,file='ph.dat',status=' unknown ' ) 
write(16. Ill ) < t ime < i > . bl < i ) . i = 1 , n> 

111 format(lx.fl2.5.8x.fl2.5> 

write (16. 211) ( t ime ( i ) . b2 ( i ) . i =1 , n> 

211 format(lx.fl2.5.8x,fl2.5> 

close(unit=16,status=’keep') 

stop 

end 

function convert ( x , i p) 

common / calib / v( 40 , 2 ) , h (40 , 2 ) , nl , n2 

i f ( ip. eq. 1 > n=nl 
i f ( ip. eq. 2) n=n2 
do 1 0 i =1 , n-1 

1 f ( X . 1 t . v( i . i p> . and . X . ge . V ( i +1 . i p> > go to 100 
10 continue 

print *. 'interpolation failure' 
print 'probe number '.ip 
print *. X 
stop 

100 grad=(h( i + l . i p) -h( i . ip) )/ ( v( i +1 . i p> -v ( i . i p> ) 

convert=h( i . i p) tgrad* (x-v(i.ip)) 


return 

end 



PROGRAMME TO CALCULATE WAVENUMBER & CELERITY FOR RANDOM WAVE 

CALCULATION OF POWER SPECTRUM OF SIGNALS USING FFT 
CALCULATION OF PHASE DIFFERENCE BETWEEN TWO SIGNALS. 

dimension al(1100>.a2<1100).b(1100) .powerl (1100) ,power2( 1 100 
dimension al 1 ( 1 1 00 ) . a22 ( 1 1 00 > 
complex xl ( 1 1 00 ) , x2< 1 100 ) 

dimension pl 2 ( 1 1 00 > . t ime < 1 1 00 ) . pcr( 1 1 00 > . pi 1 ( 1 1 00 ) . p22 < 1 1 00 ) 
dimension ak<1100>,cel<1100). phased ( 1100) 
dimension w(1100> 


pi=5.1415 
nu = l 0 
delx= .10 


n=2**nu 
pi=>. 141^926 

open(unit=8,fi le=*ph.dat'.status='oId’ > 

read(8,*> < t i me < i > . al 1 ( i > , i =1 . n) 

read<8,*) < t ime ( i ) .a22 ( i ) . i =1 , n> 

sea 1 e = l .0 

do 20 1 = 1 . n 

al< 1 >=all( 1 >*scale 

a2< 1 )=a22< i >*scale 

cont inue 

ciose(unit=8,status='keep’ > 


del=time<2>-time(l> 
tt=de i*(n-l > 

call f f t (al . b , xl , nu , de 1 ) 
call f f t (a2 , b , x2 , nu , de 1 ) 

CALCULATION OF COHERENCE 

do 100 1 = 1 .n/2 
rl = rea 1 < xl ( i > ) 
sl = imag<xl ( i ) > 
r2 = reaT < x2 < i ) ) 
s2= imag<x2 ( i > > 
pl 1 < 1 ) =atan (si /rl ) 

p22( 1 )=atan(s2/r2) 
if ( rl . 1 t . 0 . and . s 1 . gt . 0 ) 
if ( rl . 1 t . 0 . and . s 1 . It . 0 ) 
if ( rl . qt . 0 . and . si . 1 t . 0 ) 
if ( r2 . Tt . 0 . and . s2 . gt . 0 > 
if ( r2 . 1 t . 0 . and . s2 . I t . 0 > 
if ( r2 . gt . 0 . and . s2 . 1 t . 0 ) 


tl=rl*r2+sl*s2 

t2=sl*r2-rl*s2 


FUNCTION AND PHASE ANGLE. 


pl 1 ( 1 ) =pl 1 ( 1 > + pi 
pl 1 ( 1 > =pl 1 ( i > + pi 
pl 1 ( 1 ) =pl 1 ( i > + 2*p 1 
p22( 1 )=p22( i >f pi 
p22 ( I ) =p22 ( 1 > + pi 
p22( 1 >=p22( i )+ 2*pi 


pcrm=sqr t (tl**2 + t2**-2> 
powerl ( 1 > = sqrt( rl*'*f-2+sl**-2 ) 
power 2 ( i >=sqrt( r2**2+s2**2 ) 
pcr( 1 > =pcrm/ ( powerl ( i )*power2 ( i ) ) 
continue > 


do 50 i =7 , n/2 , 5 

if (pl 1 ( 1 ) . gt . pl 1 ( 1 -5 ) ) go to 50 
pll ( i )=pll ( 1 )+2*pi 
go to 25 
cont 1 nue 

do 60 i=7.n/2. 5 

if (p22( 1 > .gt .p22( i -5> > go to 60 

p22( 1 >=p22( 1 )+2*pi 



65 

555 


56 


51 

59 

222 

555 


85 

86 

1150 

1100 

1150 

1190 


do 65 1=2. n/2 . 5 
pl2(i>=abs<pll( i)-p22(i )) 
cont 1 nue 

open<un it=17,file=’ax') 

wn te( 17. 555> (pi 1 ( i ) . p22 ( i ) . pi 2 < i > , i=2 .n/2, 5> 

format (5x.fl2.5.8x,fl2.5.4x,fl2.5> 

pmax=l . Oe-08 

do 56 1=2. n/2. 5 

pmaxl=amaxl (powerl ( i ) . pmax> 

1 f ( pmax . ne . pmaxl ) ii=i 
pma X = pma x 1 
cont 1 nue 

print *. 'maximum power and location', power 1 ( i i >. i i 

do 51 1=2. n/2. 5 

phased{ i ) =pl 2 ( i > 

ak ( 1 ) =phased ( i ) /de 1 x 

cel(i>=2.0*pi*b<i )/ak(i> 

cont 1 nue 

do 59 1=2. n/2. 5 
w<i)=2*pi*b(i> 
cont 1 nue 

open(unit = 7. f i le= 'output' > 
write(7.222>(w(i).ak(i). 1=2.450.5) 
format (5x.fl2.5.8x,fl2.5) 

open(unit=15.fi le=' i nput ' > 
wr ite(15,55M<w<i>.cel( i>. i=2.450.5> 
format(5x.fl2.5.8x.fl2.5> 
print*. 'coherence factor' .pcr< i i ) 

delf=l .0/tt 

open<unit=22.f i le=' power . dat ' . status = ' unknown ' > 
fnyquist=l .0/(2. 0*de 1 > 
do 85 i =1 . n 

1 f < b< i ) . ge . f nyqu i st ) go to 86 
cont 1 nue 

write (22.1150) 

f ormat ( /2x ,' frequency in Hz; amplitude cmsquare . sec ' . // ) 
write(22. 1100) 

format (5x. 'frequency' .llx. ' ampl i tude ' , 12x . ' phase ' . / / ) 
write(22.1150) (b(i). powerl ( i ) . power2( i).pl2(i>.i=l.i)) 
format (5x.el2.5,4x,el2.5,4x.el2.5,4x,el2.5) 
wn te( 22 . 1 1 90 ) fnyquist 

format(/5x. 'Nyquist frequency=’ .el2.5) 
close(unit=22.status='keep') 


stop 

end 



6 


10 

20 


25 

55 

50 


10 


subroutine f f t (a , b . x , nu . de 1 > 
dimension xdlOO) ,b(1100>,a<1100> 
compl ex X , u . w. t 
n = 2**nu 

pi=5. 141592655 
tt=(n-l >*del 


do 5 
b( ) 1 
cont 


J ] = 1 .n 

) = < j j -1 >/tt 

1 nue 


do 6 1 = 1 , n 
X ( i ) = cmp lx(a<i>,0.0> 
cont 1 nue 


Do 20 1 =1 , nu 
1 e = 2**{ nu+1 - 1 > 
lel=le/2 
u=(l . .0. > 

w= cmp ix<cos<pi/float(iel)),-sin<pi/float(lel>)> 

do 20 ] =1 , 1 si 

do 10 1 = ) , n . 1 e 
1 p = 1 + 1 el 
t=x{i>+x(ip> 
x( ip)=(x< 1 >-x( ip>)*u 
x< i > = t 
u = u*w 


nv2=n/2 
nml =n-l 
I = 1 

do 50 1 = 1 , nml 

i f < 1 . ge . ) >go to 25 

t = x( jT 

X ( 1 ) =x< 1 > 

x( 1 )=t 

k=nv2 

1 f <k . ge . J >go to 50 


go to 55 
T=j+k 

return 

end 

subroutine trap<x,n,e> 
dimension x<1100) 

fx = 0.0 

do 10 i=2,n-l 
f x = f x + x( i > 
cont 1 nue 

e=<x( 1 >+x<n)+2 . 0*fx)/ <2 . 0*(n-l ) ) 


return 

end 



PROGRAMME TO CALCULATE WAVE SPECTRUM 


c 


c 


65 

c 


75 


45 


100 


85 

86 

1150 

1100 

1150 

1190 


dimension a(5200>, power ( 5200) ,auto<5200) .aiag( 5200) 
dimens ion b( 5200) ,a2< 5200) , 1 1 me <5200) 
complex x<5200) 
common / tr /troom 

print *, 'what are nu ndata del' 
read *. nu, ndata, del 
seal e=l . 0 

n=2**nu 
tt=de l*<n-l ) 
pi=5. 1415926 

open ( un it=21 ,f i le= ’p’ ,status='old' ) 
read(21.*) <ti me <i),a(i),i = l. ndata > 
close<unit=21 ,status='keep' ) 

remove stray mean. 

call trap<a , ndata , e) 
do 65 1 = 1 , ndata 
a< 1 ) =a( i > -e 
cont 1 nue 


filter data <Tukey- Hanning formula). 

1 1 1 =0 . 5*de 1* (ndata -1 ) 
do 75 1=1, ndata 
tl=del*< 1 - 1 ) 
wt = 1 .0 

a2 ( 1 ) =a( 1 )**2 

a( 1 )=a< 1 )*wt*scale 

cont 1 nue 


call t rap(a2 , ndata , van ance ) 

do 45 i=ndata+l,n 
a( i )=0 . 0 
continue 

call f f t (a . b , X , nu , de 1 ) 


do 100 i = 1 ,n 
r=real <x< i ) ) 
s= imag<x< i ) ) 

power ( i )=<<sqrt<r**2+s**2) >/tt) 
cont i nue 

print •*■, 'power spectrum calculated 
delf=l .0/tt 

open(uni t=22 . f i le=’pp.dat’ , status =’ unknown ' ) 

f nyqu ist=1.0/<2. 0*de 1 ) 

do 85 1 = 1 , ndata 

i f <b( i ) .ge. fnyquist) go to 86 

cont i nue 

write (22,1150) . j. j 1 

format(/2x. 'frequency in Hz ; amplitude cmsq/sec ,//> 

wr i te ( 22 , 1 1 00 ) , ,, , i 4 . ^ • //s 

format (5x. 'frequency' . 1 lx, 'ampl i tude , //) 
write<22,1150) (b(i) ,power( i>, i = l . i j > 
format (5x,el2.5,8x,el2.5) 
write(22,1190) fnyquist 

format</5x, 'Nyquist frequency= ,el2.5) 
close(unit=22,status='keep' ) 


stop 

end 



7C 


5 

6 


10 

20 

25 

55 

50 


subroutine f f t (a , b . x , nu , de 1 > 
dimension x(5200),b(5200>,a(5200) 
c omp lex X , u , w , t 
n=2**nu 

pi = 5. 141592655 
tt=(n-l )*del 

do 5 j j =1 , n 
b<);> = (jj-l)/tt 
cont 1 nue 

do 6 1 =1 , n 

x< 1 >=cmplx<a( 1 ) , 0 . 0 ) 
cont inue 

do 20 1 =1 . nu 
le=2**(nu+l-l > 
lel=le/2 
u=<l . .0. > 

w=cmpix<cos<pi/float< lei ) ) . -sin(pi/float< lei) > ) 

do 20 )=1 . lei 

do 10 1 = ) . n , 1 e 
1 p = 1 + 1 el 
t=x< 1 ) +x( ip> 
x( ip>=<x( 1 >-x< ip) >*u 
x< 1 ) =t 
u = u*w 

nv2=n/2 
nml =n-l 
1 = 1 

do 50 1 = 1 , nml 
I f ( 1 . ge . J igo to 25 

t = x ( J T 

x( j ) =x ( i > 
x( i )=t 
k = nv2 

X f ( k . ge . j >go to 50 
) = j -k 
k = k/2 
go to 55 

j = j+k 

return 

end 

subroutine trap(x.n.e> 
dimension x(5200) 

fx=0. 0 

do 10 i=2,n-l 
f x = f x+x( 1 ) 
cont 1 nue 

e=(x( 1 >+x(n)+2 . 0*fx)/(2 . 0»Cn-l ) ) 


return 

end 



400 


PROGRAMME TO CALCULATE ENERGY FLUX 

dimension cel(1000).ak<1000>,cg<1000).w(1000).s<10000> 

g=9.81 
df = .0156 
rho=l 000 
n = 86 


print*, "d=?” 
read*, d 

open <uni t = 55 , f 1 le = ’put ' ) 
read( 5?,*)<cei<i),ak<i>,i=l,n> 

do 1 =1 ,n 

cg< 1 >=0.5*cel< i )*< 1 +2*ak ( i ) *d/s inh(2*ak( i >*d> ) 
enddo 

open <un 1 1 = 44 , f 1 1 e = ' output! ' > 
write<44,400> (cel(i),cg(i),i=l,n> 
f o rma t<lx,fl2.5,8x,fl2.>) 

open(unit=55.fi le=*pp.dat’) 
read(55,*)(w<i > ,s< i) , i=l ,n) 
sum=0 . 0 

do 1 = 1,86 

sum=sum+s( i )*cg( i > 
enddo 

p r 1 n t * , " s um= " , s um 

ef lux = rho*df *g*sum 
print*, ”eflux=" ,ef lux 


stop 

end 



c 


PROGRAMME TO CALCULATE SPECTRUM WIDTH PARAMETER 


dimension f < 1 000 > , s< 1 000 > 

print*. ”df = ''.df 

read* , df 

pi=?. 1415 

dw=2*p i *df 

open ( un 1 t = 7 , f 1 1 e= ' 1 nput ' > 
read(7.*)(f ( i ) , s( i ) . i =1 . 512 > 

suml =0.0 
sum2=0 . 0 
sum?=0 . 0 

do 1 = 1,512 

suml = suml +s < i > 

sum2= sum2+ ( s < 1 )*( f ( i )**2 > ) 

sum5= sum5+ ( s < 1 > *( f ( 1 )**4 > ) 

enddo 

print*. "suml =" , suml 
print*. ”sum2=” , sum2 
print*. "sum5=" .sum? 

s 1 gmasquare = ( 1 - < sum2/ (suml *sum5*df ) > > 
pr 1 nt* , " s i gmasquare" . s i gmasquare 

stop 

end 



c 


PROGRAMME TO GENERATE REGULAR WAVE 


j^^include <default.h> 
/^include <defaul t . ttf> 
^include <macros.h> 


run "1 hz wave 0.05m" with 
OneFreq 0<1 .0,0.05) 

sample ""Vlaveguage 51 Jan 1994" on 0 
end 


c 


regular wave parameters frequency=l Hz 

ampl i tude=5 cm 



c 


PROGRAMME TO GENERATE RANDOM WAVE 


^include <default.h> 

^include <de fau 1 t . t t f > 

(^include <macros.h> 

wave wl = j onswapC 0 . 7? , 0 . 0081 , 1 . 0 , 0 . 07 , 0 . 09 ) 
wave w2 = random<wl,l) 

run "Jonswap spectrum" with 
wavemaker w2 with TTF on 0 
sample "Waveguage <P1> 01 Mar 1994" on 0 
sample "Wavegauge (P2) 01 Mar 1994" on 0 
end 


c 


jonswap parameters 


peak frequncy 

=.75 Hz 

alpha 

=.0081 

gamma 

= 1.0 

s i gma a 

= . 07 

s i gma b 

= .09 



